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Abstract. We address the question of which quantum states can be inter-converted 
under the action of a time-dependent Hamiltonian. In particular, we consider the 
problem as applied to mixed states, and investigate the difference between pure 
and mixed-state controllability introduced in previous work. We provide a complete 
characterization of the eigenvalue spectrum for which the state is controllable under 
the action of the symplcctic group. We also address the problem of which states can 
be prepared if the dynamical Lie group is not sufficiently large to allow the system to 
be controllable. 



Criteria for reachability of quantum states 



2 



1. Introduction 

The subject of control of quantum systems has been a fruitful area of investigation 
lately. The growing interest in the subject can be attributed both to theoretical 
and experimental breakthroughs that have made control of quantum phenomena an 
increasingly realistic objective, as well as the prospect of many exciting new applications 
such as quantum computers ^ or quantum chemistry 0, which attracts researchers 
from various fields. 

Among the theoretical problems that have received considerable attention lately 
is the issue of controllability of quantum systems. Various aspects such as the 
controllability of quantum systems with continuous spectra 0, wavefunction 
controllability for bilinear quantum systems P, ||, 0, controllability of distributed 
systems controllability of molecular systems 0, controllability of spin systems |[I0[] , 



controllability of quantum evolution in NMR spectroscopy [TT|, and controllability of 



quantum systems on compact Lie groups O, 13, 14, 15] have been addressed, and related 



problems such as the dynamical realizability of kinematical bounds on the optimization 
of observables |1^, [1^, and the relation between controllability and universality of 



quantum gates |T8[ , as well as the information-theoretic limits of control have been 
studied. 

In this process, various notions of controllability have been introduced. Recent work 
on controllability of quantum systems on compact Lie groups has finally shown that the 
degree of controllability of a quantum system depends on its dynamical Lie group, 
and that many different notions of controllability are in fact equivalent. In particular, 
it has been proved that quantum systems evolving on a compact Lie group, such as 
closed quantum systems with a discrete energy spectrum, are either density matrix 
/ operator controllable, pure-state / wavefunction controllable, or not controllable 
Il20| , pl| . For density matrix, operator or completely controllable quantum systems, every 
kinematically admissible target state or operator can be dynamically realized, and the 
kinematical bounds on the expectation values (ensemble averages) of observables are 
always dynamically attainable [17]. Fortunately, many quantum systems have been 



shown to be completely controllable ]13, 14, pT] 



Nevertheless, there are quantum systems that are either only pure-state controllable 
or not controllable at all. For instance, it has been shown that the dynamical Lie group 
of certain atomic systems with degenerate energy levels is the (unitary) symplectic 
group, which corresponds to pure-state controllability ]|21|. Other systems with certain 
symmetries may be either pure-state controllable or non-controllable depending on the 
symmetry. For instance, given a system with equally spaced energy levels and uniform 
dipole moments for transitions between adjacent levels, the dynamical Lie group is the 
symplectic group if the dimension of its Hilbert space A'^ is even, but it is the orthogonal 
group if is odd [^. For these systems, the question of dynamical reachability of target 



states, which is important in many applications, remains. In this paper, we address this 
problem by studying the action of the djTiamical Lie group of pure-state-only and non- 
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controllable quantum systems on the kinematical equivalence classes of states. Explicit 
criteria for dynamical reachability of states are derived for systems whose dynamical Lie 
group is the (unitary) symplectic group or the orthogonal group. 

2. Quantum states and kinematical/dynamical equivalence classes 

We consider a quantum system whose state is represented by a density matrix acting 
on a Hilbert space Ti. of dimension N. A density matrix always has a discrete spectrum 
with non- negative eigenvalues Wn that sum to one, J2n — and a spectral resolution 
of the form 

TV 

p=J2^n\'^n){'^n\, (l) 
n=l 

where are the eigenstates of p. The for 1 < n < N are elements of the Hilbert 
space Ti. and can always be chosen so as to form a complete orthonormal set for Ti.. The 
{'^n\ are the corresponding dual states defined by 

I = 5mn Vm, n. (2) 

Conservation laws such as conservation of energy and probability require the time 
evolution of any (closed) quantum system to be unitary. Thus, given a Hilbert space 
vector l^o)) its time evolution is determined by |^(t)) = C/(i)|^o) where U{t) a unitary 
operator for all t and U (0) = /. Hence, a density matrix po must evolve according to 

p{t) = U{t)poU{t)\ (3) 

where U{t) is unitary for all times. This constraint of unitary evolution induces 
kinematical restrictions on the set of target states that are physically admissible from 
any given initial state. 

Definition 1 Two quantum states represented by density matrices po and pi are 
kinematically equivalent if there exists a unitary operator U such that pi — UpoW. 

Thus, the constraint of unitary evolution partitions the set of density matrices on 
Ti into (infinitely many) kinematical equivalence classes. It is well known that two 
density matrices po and pi are kinematically equivalent if and only if they have 
the same eigenvalues. The kinematical equivalence classes are therefore determined 
by the eigenvalues of p. Furthermore, we introduce the following classification of 
density matrices according to their eigenvalues, which we shall relate to the degree 
of controllability of the system. 

Definition 2 (Classification of density matrices) Every density matrix is of one 
of the following types. 

(i) Completely random ensembles: Density matrices whose spectrum consists of a 
single eigenvalue Wi — jr that occurs with multiplicity N. 
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(ii) Pure-state-like ensembles: Density matrices whose spectrum consists of two distinct 
eigenvalues, one of which occurs with multiplicity one and the other with multiplicity 
N -I. 

(Hi) General ensembles: Density matrices whose spectrum consists of at least two distinct 
eigenvalues, at least one of which occurs with multiplicity Ni where 2 < Ni < N — 2; 
or density matrices whose spectrum consists of N distinct eigenvalues (N >2). 

Note that type (0) (pure-state-like ensembles) includes density matrices 
representing pure states such as p = diag(l, 0, 0, 0) but not every density matrix in 
this class represents a pure state. For instance, p = diag(0.7, 0.1, 0.1, 0.1) is of type (P) 
but does not represent a pure state. 

Given a specific quantum system with a control-dependent Hamiltonian of the form 

M 

H[fiit), . . . , /mW] = Ho+J2 Umm, (4) 

m=l 

where the fm, 1 < m < M, are (independent) bounded measurable control functions, the 
question arises which states are dynamically reachable from a given initial state. Clearly, 
the set of potentially dynamically reachable states is restricted to states within the same 
kinematical equivalence class as the initial state. However, not every kinematically 
admissible target state is necessarily dynamically reachable. Since the time-evolution 
operator U (t) has to satisfy the Schrodinger equation 

injU{t) = H[f,{t),...jM{m{t), (5) 

where H is the Hamiltonian defined above, only unitary operators of the form 

U{t) = exp+ l^-^H [/i(t), . . . , /M(t)]} , (6) 

where exp_,_ denotes the time-ordered exponential, qualify as evolution operators. Using, 
for instance, the Magnus expansion of the time-ordered exponential, it can be seen that 
only unitary operators of the form exp(x), where x is an element in the dynamical Lie 
algebra L generated by the skew- Hermit ian operators iHq, . . . ,iHM, are dynamically 
realizable. These operators form the dynamical Lie group S of the system. 

Definition 3 Two kinematically equivalent states po and pi are dynamically equivalent 
if there exists a unitary operator U in the dynamical Lie group S such that pi = U p^U'^ . 

This dynamical equivalence relation subdivides the kinematical equivalence classes. 

In the following, we shall be particularly concerned with the unitary group U{N), 
the special unitary group SU{N), the (unitary) symplectic group Sp^^) the 
(unitary) orthogonal group SO{N). As usual, the unitary group U{N) is the compact 
Lie group consisting of all regular N x N matrices U that satisfy WU = UW = I. The 
special unitary group SU {N) is the subgroup of U{N) consisting of all unitary matrices 
U G U{N) whose determinant is +1. For our purposes in this paper, we define the 
symplectic group and the special orthogonal group as follows. 
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Definition 4 The (unitary) symplectic group Sp{l) is the subgroup of SU{2£) 
consisting of all unitary operators of dimension 2£ that satisfy JU = J for 

where is the identity matrix of dimension i. 

Definition 5 The (unitary) special orthogonal group SO{N) is the subgroup of SU{N) 
consisting of all unitary operators of dimension N that satisfy U^JU = J for 

/ T \ / 1 \ 

J= r N = 2i, J = 



0/, 
V /, y 



N = 2i+1. (8) 



3. Dynamical Lie group action on the kinematical equivalence classes 

The set of quantum states that is dynamically accessible from a given initial state po 
depends on the action of the dynamical Lie group S on the kinematical equivalence 
classes of density operators. 

Definition 6 The dynamical Lie group S of a quantum system is said to act transitively 
on a kinematical equivalence class C of density matrices if any two states in C are 
dynamically equivalent . 

Since the equivalence class of completely random ensembles [type (0) above] consists 
only of a single state p = j^In, it follows immediately that every group acts transitively 
on this equivalence class. 

Any dynamical Lie group 5" that does not act transitively on the kinematical 
equivalence class of pure states, acts transitively only on the trivial kinematical 
equivalence class of completely random ensembles. Furthermore, from classical results 
by Montgomery and Samelson [^, it follows that U{N), SU{N), Sp{^N) and 
Sp{^N) X U{1) are the only dynamical Lie groups (up to isomorphism) that act 
transitively on the equivalence class of pure states. Therefore, any dynamical Lie 
group S that is not isomorphic to either U{N), SU{N), Sp^^N) or Sp^N) x U{1) 
acts transitively only on type (|) states, i.e., completely random ensembles. U{N) and 
SU{N) clearly act transitively on every kinematical equivalence class of states, which 
leaves only Sp{^N) and Sp{^N) x U{1), whose action on the kinematical equivalence 
classes of states we shall now address. 

We begin by showing that transitive action of Sp{^N) on pure states implies 
transitive action on all equivalence classes of type We shall prove this result for the 
standard representation of Sp{^N) as defined above. To see that this is sufficient, note 
that lemma 4.2 in [EO] shows that whenever the dynamical Lie algebra of a quantum 



system of the type considered in this paper is isomorphic to sp{^N), then it is conjugate 
to sp{^N) via an element in U{N). Thus, if the dynamical Lie group S of the system is 
of type Sp{^N) then it is not only isomorphic to the standard representation of Sp{^N), 
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but there exists a unitary transformation (basis change) B that maps any unitary 
operator m U G 5 to a unitary operator U = BUB'^ in the standard representation 
of Sp{\N), i.e., S is unitarily equivalent to the standard representation of Sp{\N). 

Note that theorem 6 in |]20[ gives a general condition for transitive action of a 
dynamical Lie group S C U{N) on a kinematical equivalence class of states represented 
by a density matrix p: the action is transitive if and only if 

dim U{N) - dim 5 = dim Cp - dim{Cp n S) , (9) 

where Cp is the centralizer of p and Cp fl is the intersection of the centralizer with S. 
However, since determination of the dimension of Cp, and especially Cp fl S, tends to be 
very difficult in practice (see [Appendix B| for an example) we shall not use this result 



but pursue an alternative approach instead. 

Lemma 1 Sp{^N) acts transitively on all kinematical equivalence classes of density 
matrices whose eigenvalues satisfy Wi ^ W2 = w^, = ■ ■ ■ = w^- 

Proof: Any p with eigenvalues Wi W2 = = . . . = wn can be written as 

p = Wi\<i/){^\+UJ2P{\'^)^), 

where P(|\E')-'") is the projector onto the orthogonal complement of the subspace spanned 
by Hence, any pair of kinematically equivalent states of this type is of the form 

P0 = ^l|^(0))(^(0)|+^2P(|^(0))±) 
= + ^«2P(|^(1))^). 

Since Sp{^N) acts transitively on the equivalence class of pure states, there exists a 
unitary operator [/ G Sp{^N) such that Ul"^^^^) = |\E'(^)). Furthermore, U automatically 
maps the orthogonal complement of 1^^^°^) onto the orthogonal complement of 
since it is unitary and thus we have 



P 



(1) 



Hence, Sp{^N) acts transitively on all equivalence classes of density matrices whose 
eigenvalues satisfy wi ^ W2 = = . . . = wn- ■ 

However, the action of Sp{^N) on the class of pure states is not two-point transitive 
as the following example shows. 

Example 1: Let N = 2i and a and b be two unit vectors in C^. Since = 2i, we 
can partition the vectors as follows 





where aj, bj for j = 1,2 are vectors in C . Since Sp{i) acts transitively on the unit 



sphere in (D^ it follows that there exists a U E Sp{i) such that Ua = b. However, since 
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any unitary operator in Sp{£) satisfies U^JU = J with J as in (|^, we have U = .PU*J 
and thus J'^U* Ja = b or equivalently UJa* = Jb*. Noting that 

So \ .r. f -b% 



Jo* =[ , Jb* 



2 



b\ 



it thus follows that U maps c = J a* onto d = Jb*. Therefore, given two (orthogonal) 
unit vectors of the form a and c, it is not possible to find a unitary transformation in 
Sp{i) that maps these two vectors onto two arbitrary (orthogonal) unit vectors. Rather, 
once we have chosen the image of a, the image of c is fixed. ■ 

This lack of two-point transitivity has serious implications for the action of Sp{^N), in 
particular it implies non-transitive action on all kinematical equivalence classes of type 



Lemma 2 Sp{^N) does not act transitively on kinematical equivalence classes of 
density matrices with at least three distinct eigenvalues, two of which having multiplicity 
one. 

Proof: Any two kinematically equivalent density matrices can be written as 

N N 

P0 = ^m„|^'„)(^'„|, pi = Y.Wn\'^n){'^n\- 

n=l n=l 

Since there are at least three distinct eigenvalues and two of them have multiplicity one, 
we may assume Wi ^ Wn for all n 7^ 1 and W2 7^ Wn for all n 7^ 2. Thus, |\E'„) and 
for n = 1, 2 are unique up to phase factors and any U such that pi = UpoW must map 
onto (modulo phase factors) for n = 1,2, i.e.. 

However, suppose |\E'i) = a, |\&2) = c and |$i) = 6 but |$2) 7^ e^'^rf, where a, b, c and 
d are as defined in example |I|. This example then shows that it is impossible to find a 
U G Sp{^N) that simultaneously maps |\E'i) onto and |\E'2) onto |$2)- Therefore, 
there does not exist a unitary operator in Sp{^N) such that pi = U p^U^ . ■ 



Lemma 3 Sp{\N) does not act transitively on equivalence classes of density matrices 
that have at least one non-zero eigenvalue that occurs with multiplicity greater than one 
but less than N — 1. 

Proof: Suppose wi has multiplicity A^i where 2 < A^i < - 2. If Sp{\N) acts 
transitively on the selected equivalence class of states then we must be able to map the 
A^i-dimensional eigenspace E^^\wi) for po onto the corresponding eigenspace E'^^\wi) 
for pi by a unitary operator in Sp{\N) for any po and pi in the same equivalence 
class. However, it is easy to see that this is not always possible. Suppose E^'^\wi) 
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contains a pair of vectors of the form defined above and E'^^\wi) contains a 

vector 6 but the related vector d is in the orthogonal complement of E^^\wi). Then it 
is impossible to map [wi] onto E^^\wi) by a f/ G Sp{-^N). Since the orthogonal 
complement of E^^'^^wi) has at least dimension two, we can always choose E^^\wi) such 
that d G E^^\wi)-^. Hence, Sp{^N) does not act transitively on the selected equivalence 
class of states. ■ 

Given any two mixed states po and pi related by pi = UpoW for some U G 
Sp{^N) X f/(l), we can find a [/ G Sp{^N) such that pi = UpoW. For instance, if 
det U = e*°, setting f/ = e'^^^^U produces an operator with det{U) = 1 that obviously 
satisfies 

UpoU^ = UpoU^ = Pi. 

Thus, Sp{^N) X U{1) acts transitively on a kinematical equivalence class C of density 
matrices if and only if Sp{^N) does. Combining this observation with the previous 
lemmas yields the following theorem. 

Theorem 1 

• U{N) and SU{N) act transitively on all kinematical equivalence classes. 

• Sp{^N) and Sp{^N) x U{1) act transitively on all kinematical equivalence classes 
of density matrices of type (§j or ^ and only those. 

• Any other dynamical Lie group acts transitively only on the trivial kinematical 
equivalence class of completely random ensembles. 

4. Criteria for reachability of target states 

Having established that the action of the dynamical Lie groups Sp{^N) and Sp{^N) x 
f/(l) is not transitive on any kinematical equivalence class of density matrices of type 
([!]), and that all other dynamical Lie groups except U{N) and SU{N) act transitively 
only on the trivial kinematical equivalence class of completely random ensembles, the 
question of identifying states that are kinematically but not dynamically equivalent 
arises. 

Since dynamical Lie groups can be very complicated, it would be unrealistic to 
expect that simple criteria for dynamical equivalence of states can be derived for 
arbitrary dynamical Lie groups. However, for certain types of dynamical Lie groups 
of special interest, such as Sp^N) [or Sp{^N) x U{1)] and SO{N) [or SO{N) x [/(I)], 
this is possible, as will be shown in the following. 

4.1. Systems with dynamical Lie group Sp{^N) or Sp{jN) x f/(l) 

To address the problem of finding criteria for dynamical equivalence of states for systems 
whose dynamical Lie group S is isomorphic (unitarily equivalent) to Sp{^N), we recall 
that any unitary operator U G Sp{i) satisfies U^JU = J for J as defined in (^. Thus, 
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any dynamical evolution operator U for a system of dimension N = 2i with dynamical 
Lie group of type Sp{i) must satisfy 

U^JU = J (10) 

for a matrix J, which is unitarily equivalent to Therefore, we must have 

U = J^U*J, = J^U^J. 

Two kinematically equivalent states po and pi are thus dynamically equivalent if and 
only if there exists a unitary operator U such that 

= UpoU^ and pi = J^U*JpoJ^U^J, 

or equivalently, 

pi = UpoU^ and {JpiJ^y = U ( Jpo J^)* Ul (11) 



pi PO 



Example 2: Let = 4 and 5" = 5*^(2) with J = J as in (0). 

(i) Then po = diag(a, a, 6, 6) {0 < a,b < ^, a + b = ^) and pi = diag(a, fo, 6, a) are 
dynamically equivalent since there exists a unitary operator U such that pi = U p^U^ 
and any such U clearly maps po = diag(6, 6, a, a) to pi = diag(6, a, a, 6). 

(ii) Po and p2 = diag(a, 6, a, 6), on the other hand, are not dynamically equivalent 
(unless b = a) since p2 = P2 but po ^ po and there cannot be a unitary operator 
such that pi = U pqU'^ = UpoW if po 7^ po- 

This shows that S = Sp{2) divides any kinematical equivalence class of states with two 
distinct eigenvalues of multiplicity i = 2 into at least two disjoint subsets of dynamically 
equivalent states. ■ 

Sometimes the condition U^JU = J can also be used directly to show that two states 
are not dynamically equivalent. 

Example 3: Consider again = 4 and S = Sp{2) with J = J as in (0) as well 
as the initial state po = diag(a, b, c, d) where < a,b,c,d < 1, a + b + c + d= 1 and 
a, b, c, d mutually different. We can conclude that the state pi = diag(6, a, c, d) is not 
dynamically equivalent to po since we would require a unitary operator of the form 



U 



( e^"^! \ 

e*<^2 



V e"*"^ ) 



which does not satisfy U'^ JU = J. ■ 

X See Appendix A| for details about how to determine J. 
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Another way of showing that two (kinematically equivalent) density matrices are not 
dynamically equivalent is to prove that (|Tl]) cannot have a solution by showing that the 
related linear system 

Pif/-?7po = 0, p,U-Upo = (12) 

does not have a solution. To verify this, we note that the linear system above can be 
rewritten in the form AU = where ^ is a matrix with 2N'^ rows and N"^ columns and 
f/ is a column vector of length A^^. If the null space of A is empty then there is no U 
such that AU = and hence there is no x matrix U that satisfies ([I^). However, 
note that if the linear system above does have a solution, this does not imply that the 
states in question are dynamically equivalent since the solution to the linear equation 
is in general not unitary. 



4-2. Systems with dynamical Lie group SO{N) or SO{N) x U{1) 



From the previous discussion, we know that SO{N) does not act transitively on any 
kinematical equivalence class other than the trivial one. However, we can establish 
criteria for dynamical equivalence of states similar to those for Sp{^N) by noting that 
any unitary operator U G SO{N) must satisfy U^JU = J for J as in (|^). Therefore, two 
kinematically equivalent states Pq and pi are dynamically equivalent under the action 
of a dynamical Lie group S which is unitarily equivalent to SO{N), if there exists a 
unitary operator U such that 

UpoU^ and {JpiJ^y = U {JpoJ^Y Ul (13) 



Pi 



PO 



with J unitarily equivalent to (||), and determined as described in [Appendix A . 
Example 4: Consider a system with N = 5 and Hamiltonian H = Hq + f{t)Hi where 



Ha 



/ 


-2 











\ 
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2 j 
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1 





It can be verified using the algorithm described in that the Lie algebra of this 
system has dimension 10, which is equal to the dimension of so(5). Using the technique 
described in [Appendix A| , we find that both of the generators iHq and iHi of the Lie 
algebra satisfy J + Jx = for 



J 



( ° 











1 \ 











-1 











1 











-1 











1 1 











J 
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which is unitarily equivalent to the standard J for so(5). We can thus conclude that 
its dynamical Lie algebra is so(5) and its dynamical Lie group is 5*0(5). Furthermore, 
note that the two pure states 



/ 1 
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^ 0.5 











0.5 y 



are not dynamically equivalent since {JpiP)* = pi but {Jp^P)* ^ po and it is thus 
impossible to find a unitary transformation such that UpoW = pi = U{Jpo,r)*W . ■ 



5. Conclusion 

The question of dynamical equivalence of kinematically equivalent quantum states has 
been been addressed by studying the action of the dynamical Lie group of the system on 
the kinematical equivalence classes. For systems whose dynamical Lie group is unitarily 
equivalent to either Sp{^N) or SO{N), explicit criteria for dynamical reachability / 
equivalence of states have been given, and their application illustrated with several 
examples. 

Furthermore, we have provided a classification of density matrices according to their 
eigenvalues, which divides mixed quantum states into three main types: (i) completely 
random ensembles, (ii) pure-state-like ensembles, and (iii) general ensembles. We have 
also proved that the dynamical Lie group Sp{^N) acts transitively on all equivalence 
classes of quantum states of type (D and (0), but only those. 

Although it is known that a pure-state controllable system whose dynamical Lie 
group S is isomorphic to Sp{^N) is not density matrix controllable in general |^, this 
result shows that there are more than just a few examples of kinematically equivalent 
density matrices that are not dynamically reachable from one another in this case. In 
fact, the action of S is not transitive on almost all kinematical equivalence classes. This 
is in marked contrast to the action of 5* for a density matrix controllable system, which 
is transitive on all kinematical equivalence classes, as well as the action of S for a non- 
controllable system, which is transitive only on the trivial kinematical equivalence class 
of completely random ensembles. 
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Appendix A. Finding J for dynamical Lie groups of type Sp{\N) or SO{N) 

For the results of the previous sections to be truly useful, we must also address the 
question of how to determine the J matrix of a given system. To this end, note that the 
elements of the dynamical Lie algebra L associated with the dynamical groups Sp{^N) 
and SO{N) must satisfy a relation similar to the one satisfied by the elements of the 
group, namely any x & L must satisfy 

x^J+Jx = 0, (A.l) 

where J is the same as for the related group. Thus, given a system with total 
Hamiltonian @, this implies in particular that the generators iif^ of the dynamical 
Lie algebra must satisfy ( |A.1D . 

Equation (|A.1| ) can be written as a system of linear equations of the form 



LmJ = 0, < m < M, 

where L^ is a square matrix of dimension A^^ determined by the generators iHm and J 
is a column vector of length A^^. The solutions J of the above matrix equation can be 
found by computing the null space of the operator 

( \ 

V La/ / 

If the dynamical Lie group is of type Sp{^N) or SO{N) then the nullspace contains a 
single element J, which can be rearranged into a square matrix whose eigenvalues agree 
with whose of the standard J for the group defined above. That is, concretely, 

• if N = 2i and J has two distinct eigenvalues +i and — i, both of which occur with 
multiplicity i then the dynamical Lie group is Sp{i)] 

• if N = 2£ and J has two distinct eigenvalues +1 and —1, both of which occur with 
multiplicity i then the dynamical Lie group is S0{2i); 

• if N = 2i + 1 and J has two distinct eigenvalues +1 and —1, occurring with 
multiplicity i + 1 and i, respectively, then the dynamical Lie group is S0{2i + 1); 

Hence, the algorithm not only determines J but it also allows us to decide whether the 
dynamical Lie group is of type Sp{^N) or SO{N). 

Note that the dynamical Lie group S can only be Sp{^N) or SO{N) if all the 
partial Hamiltonians Hm of the system have zero trace. However, if any of the partial 
Hamiltonians Hm has non-zero trace then the dynamical Lie group of the system can 
still be Sp{^N) x U{1) or SO{N) x f/(l). To deal with this situation, we note that 
S ~ Spl^N) X [/(I) or 5 ~ SO{N) x f/(l) is possible only if the generators 

Xm = iHm-j^ Tt{HJIn, < m < M (A.2) 
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of the related trace-zero Lie algebra L' satisfy ( [A.l| ) for < m < M and we can thus 
proceed as above to determine J . 



Appendix B. Comparison of Theorem 1 with Theorem 6 in |20 



To demonstrate the difficulty in using theorem 6 in |2y] to verify whether the dynamical 
Lie group of a system acts transitively on an equivalence class of density operators, 
we shall consider a simple example. 

Assume the dynamical Lie group of the system is Sp{2) C U{4). According to 
theorem 1 above, Sp{2) does not act transitively on the kinematical equivalence class 
represented by p = diag(a, a, b, b) with < a, < | and a + b = \ since p is of type (jiil). 



To show that the action is not transitive using theorem 6 in pO[, we note 
first that dim[/(4) = 16 and dimS'p(2) = 10. Thus, the left hand side in (||) is 
dim U{A) - dim Sp{2) = 6. 

To compute the right hand side, we need to determine the centralizer Cp of p. 
Noting that 

\ 




hi 



2 



where I2 is the indentity matrix in dimension 2, we see that p commutes with every 
unitary matrix of the form 

' \ 



Uc 



A 



2 



where Ai and A2 are arbitrary unitary matrices in U{2). Thus, the centralizer of p is 
U{2) X U{2) and its dimension is 4 + 4 = 8. 

To compute the intersection of Cp with S = Sp{2), we recall that any matrix in 
Sp{2) must preserve J as defined in (^. Concretely, this means UqJUc = J, i.e.. 

Thus, we must have A'^A2 = l2- Noting that Ai and A2 are unitary, this is only possible 
if Ai = A2, i.e., if Ai is the complex conjugate of A2, since (742)^^2 = A2A2 = I. Hence, 
the intersection of the centralizer Cp with S = Sp{2) is U{2), which has dimension 4. 
Hence, dim Cp — dim(Cp fl Sp{2)) = 8 — 4 = 47^ 6, i.e., the left and right hand side in 
(y) are not equal. Thus we have shown using theorem 6 that the action of Sp{2) on the 
kinematical equivalence class of p is not transitive. 

References 

[1] A. Ekert, P. Hayden, H. Inamori, and D. K. L. Oi, Int. J. Mod. Phys. A 16, 3335 (2001). 
[2] H. Rabitz, R. de Vivie-Riedle, M. Motzkus, and K. Kompa, Science 288, 824 (2000). 
[3] T. J. Tarn, J. W. Clark, and D. G. Lucarelli, in 39th IEEE CDC Proceedings (Causal Productions, 
Adelaide, Australia, 2000), pp. 2803-2809. 



Criteria for reachability of quantum states 



14 



[4] G. M. Huang, T. J. Tarn, and J. W. Clark, J. Math. Phys. 24, 2608 (1983). 
[5] G. Turinici, Comptes Rendus Acad. Sci. Ser. I-Math. 330, 327 (2000). 
[6] G. Turinici and H. Rabitz, Chem. Phys. 267, 1 (2001). 

[7] G. Turinici, in 39th IEEE CDC Proceedings (Causal Productions, Adelaide, Austraha, 2000), pp. 

1364-1369. 

[8] J. L. Lions, Proc. Natl. Acad. Sci. U. S. A. 94, 4828 (1997). 
[9] V. Ramakrishna et al, Phys. Rev. A 51, 960 (1995). 
[10] F. Albertini and D. D'Alessandro, quant-ph/0106115 (2001). 

[11] T. S. Untidt, S. J. Glaser, C. Gricsingcr, and N. C. Nielsen, Mol. Phys. 96, 1739 (1999). 
[12] D. D'Alessandro, in 39th IEEE CDC Proceedings (Causal Productions, Adelaide, Austraha, 2000), 
pp. 1086-1091. 

[13] H. Fu, S. G. Schirmcr, and A. I. Solomon, J. Phys. A 34, 1679 (2001). 
[14] S. G. Schirmer, H. Fu, and A. I. Solomon, Phys. Rev. A 63, 063410 (2001). 
[15] C. Altafini, quant-ph/01 10147 (2001). 

[16] J. V. Leahy and S. G. Schirmcr, in 39th IEEE CDC Proceedings (Causal Productions, Adelaide, 

Australia, 2000), pp. 1358 1363. 
[17] S. G. Schirmer and J. V. Leahy, Phys. Rev. A 63, 025403 (2001). 
[18] V. Ramakrishna and H. Rabitz, Phys. Rev. A 54, 1715 (1996). 
[19] H. Touchcttc and S. Lloyd, Phys. Rev. Lett. 84, 1156 (2000). 
[20] F. Albertini and D. D'Alessandro, quant-ph/0106128 (2001). 
[21] S. G. Schirmer, J. V. Leahy, and A. L Solomon, quant-ph/0108114 (2001). 
[22] S. G. Schirmer, L C. H. PuUen, and A. L Solomon, J. Phys. A 35, 2327 (2002). 
[23] D. Montgomery and H. Samelson, Ann. Math. 44, 454 (1943). 



